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ABSTRACT 



We present 1.3-mm subarcsecond SMA observations of the prototypical Class protostar VLA1623. We report the detection of 1.3- 
mm continuum emission both from the central protostellar component VLA1623 and two additional sources, Knot-A and Knot-B, 
which have been already detected at longer wavelengths. Knot-A and Knot-B are both located along the western cavity wall opened by 
the protostellar outflow from VLA1623. Our SMA observations moreover show that these two continuum sources are associated with 
bright, high-velocity I2 CO(2-l) emission, slightly shifted downstream of the outflow propagation direction with respect to the 1.3- 
mm continuum emission peaks. The alignment of Knot-A and Knot-B along the protostellar outflow cavity, the compactness of their 
1.3-mm continuum emission and the properties of the associated CO emission suggest that these two sources trace outflow features 
due to shocks along the cavity wall, rather than protostellar objects. While it was considered as one of the best examples of a close 
protobinary system so far, the present analysis suggests that the prototypical Class 0, VLA1623, is single on the scales a > 100 AU 
probed by our SMA observations. Moreover, we present here the second robust case of compact millimeter continuum emission 
produced by interactions between the protostellar jet and the envelope of a Class protostar, which suggests a high occurrence of 
these outflow features during the embedded phase. 
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Understanding the first steps of protostar formation is a major 
unsolved gap in modern astrophysics. Observationally, the key to 
constraining protostar formation models lies in high-resolution 
studies of the youngest protostars. Class objects, originally 
discovered at (sub)millimeter wavelengths are believed to be 
such very young accreting protostars ( |Andre etaL 1993 2000). 



action of the outflowing gas with the surrounding circumstellar 
medium ( Moraghan et al. 2008 1. 

VLA 16234-2417 (hereafter VLA1623) was the first Class 



object to be identified (Andre et al. 1990 1993 i, and it is con 



Because they are observed only / < 10 3 yr after their formation, 
while most of their mass is still in the form of a dense envelope, 
Class protostars are likely to retain detailed information on the 
initial conditions and physics of the collapse phase (see review 
by Andre et al. 2000). It is, moreover, believed that the Class 
phase represents a pivotal stage in star formation during which 
the angular momentum problem (Bodenheimer 1995 ) is solved 
through ejection processes and/or the fragmentation into multi- 
ple systems. Indeed, it has been shown that Class protostars 
drive powerful, highly-collimated protostellar jets (Bontemps 
[eTaTj [T9961 |Gueth & Guilloteau| [T9991 |Nisini et al.| |2002) . 
These jets are mostly detected through the swept-up gas form- 
ing a molecular outflow, or through a trail of high speed emission 
knots such as molecular bullets (Bach ilTer et al.|[T991[ ). The ex- 
act nature of these knots is still unclear because they could result 
from internal shocks within the jet or be induced by the inter- 



sidered to be the prototype of the class. The protostaiQwas orig- 
inally identified as a strong 1.2-mm continuum source (Andre 
19931) associated with a large-scale CO outflow (Andre 
while free-free radio continuum emission, likely 



et al. 



et al. 



1990 



to be tracing the shock-ionized base of the jet, was also re- 
ported fLeous et al. 1991 1. High-resolution VLA observations 
of this object were carried out by Bontemps & Andre] ( |1997| i, 
revealing a string of three candidate HH-like objects (labeled as 
knots A, B, and C along a western direction) closely aligned 

1993) 



along the bipolar flow traced by the CO (Andre et al 



and H2 emission (|Dent et aL 1995 1. As part of a subarcsecond 



* The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 



dust continuum survey of embedded young stellar objects with 
BIMA at 2.7-mm, |Looney et al.| ( |2000| ) found two components 
separated by 1.1"(~ 160 AU) at the position of VLA1623. The 
eastern 2.7-mm BIMA component is coincident with the centre 
of the VLA 1623 envelope as traced by the larger scale OVRO 
2.7-mm map ( |Bontemps & Andre 1997 ). The 2.7-mm western 
component is coincident with the first (Knot-A in Bontemps & 
Andre|[T997| ) of the VLA cm knots. While the eastern compo- 
nent, named VLA1623A, is agreed by all to be protostellar in 
nature, and it most likely traces the central protostellar object of 



1 In the following, we refer to the system composed of the central 
protostellar object plus envelope as a protostar. 
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VLA1623, the nature of the secondary component (VLA1623B 
in the naming system by |Looney et al. 2000| previously VLA 



source labeled as Knot-A) remains unc 



Table 1. Properties of the SMA 1.3-mm continuum sources de- 
tected 



ear. Indeed, Bontemps 



& Andre ( 1997 1 proposed that it is the first in a series of knots 
in the jet from VLA 1623, while Looney et al. ( 2000| l argued 



that it is a binary component of the VLA 1623 protostar. This 
controversy over whether VLA1623 is a binary protostar has re- 
cently been discussed at the light of 7-mm VLA observations by 
Ward-Thompson et al. (2011). While they suggest that the mil- 



limeter fluxes should mainly be due to dust continuum emission 
toward both sources, they could not firmly determine the nature 
of VLA1623B/Knot-A. In the following, we adopt the first nam- 
ing system introduced by Bontemps & Andre| ( |1997| ) to refer to 
the sources in the VLA 1623 environment. 

2. SMA observations and data reduction 

We obtained 1.3-mm high angular resolution maps of the proto- 
type Class protostar VLA1623, with the Sub-Millimeter Array 
(SMA) interferometer ( |Ho et al.] |2004). Observations were car- 
ried out on 2009 July 3 under good weather conditions (r-0. 1), 
with the Very Extended configuration (longest baseline ~ 509 m) 
of the array. 

The central frequency was set to 220.22 GHz in the lower 
side band (LSB), while the upper side band was centered at 
230.3 GHz. High-resolution spectral windows (0.2 MHz chan- 
nels) were chosen to cover the ~100kms~' around the rest 
frequencies of 12 CO(2-l) at 230.538 GHz, 13 CO(2-l), and 
C ls O(2-l). Low-resolution spectral windows (3.25 MHz chan- 
nels) were used to cover the rest of the bandwidth. The contin- 
uum was built from both the USB and LSB low-resolution spec- 
tral windows, excluding a ~50kms~ 1 spectral window around 
the e-CH 3 OH(8(-l, 8)- 7( 0, 7)) rest frequency (229.8GHz), 
which resulted in an effective line-free bandwidth of 3.2 GHz 
around a central frequency -226 GHz. The baseline solution ob- 
tained for Very Extended Configuration on 2009 July 9 was ap- 
plied, and we performed absolute flux calibration using Callisto 
observations, while bandpass calibration was obtained by ob- 
serving the bright quasar 1924-292. Phase and amplitude cal- 
ibrations during the track were carried out using the quasars 
1517-243 (primary calibrator located 15° from VLA1623, flux 
of 1.4 Jy at the time of the observations) and 1625-254 (sec- 
ondary calibrator located 1° from VLA1623, flux of 0.8 Jy at 
the time of the observations). Less than 8% of the data was 
flagged during the standard reduction procedure, using the soft- 
ware package MIR dedicated to SMA data |QT||201 1| ). By com- 
paring the fluxes of the quasars in our calibrated maps with their 
tabulated fluxes in the SMA database, the final uncertainty in the 
flux scale was estimated to be less than 15%. 

The data were later imaged (inversion and cleaning) using 
the MAPPING software (part of the GILD AS0 package). The 
primary beam FWHM of SMA at 225 GHz is -50", and pri- 
mary beam attenuation has been corrected. Natural weighting 
was used for imaging the data, leading to a final full-width at 
half maximum (FWHM) synthesized beam-size of 0.7" x 0.4", 
at a position angle of 31°. The 1-sigma rms noise on the fi- 
nal continuum map is 2. 1 mJy/beam. We smoothed the original 
spectral data by binning 2 consecutive channels, so that the fi- 
nal spectral maps have a spectral resolution of 0.5 kms" 1 . The 
typical 1-sigma rms noise level in the final CO maps is ~ 80 to 
120 mJy/beam per channel. 



2 Grenoble Image and Line Data Analysis System, software provided 
and actively developed by IRAM (http://www.iram.fr/IRAMFR/ 
|GILDAS| l 



Source 
a (J2000) (1) 
6 (J2000) (1) 
FWHM (2) 
P.A. (2 > 

Fpeak (mJy/beam) 
F int (mjy) 



VLA 1623 



Knot-A 



Knot-B 



16:26:26.39 
-24:24:30.7 
1.2"xl.l" 
20° 
59±8 
175+25 



16:26:26.31 
-24:24:30.5 
0.57"x0.24" 
14° 
71+7 
93±20 



16:26:25.63 
-24:24:29.5 
0.72"x0.27" 
13° 
26±4 
29±6 



(1) The accuracy on the position of the sources, estimated from ab- 
solute positional accuracies of the calibrators observed, is ~0.03". 

(2) Deconvolved sizes and position angles of the sources (from the 
modeling of our SMA visibilities, see Fig.[2]). 



Here, we only discuss the content and analysis of the contin- 
uum and 12 CO(2-l) maps, because interferometric filtering (the 
shortest baseline is 68 m) and low signal-to-noise ratio in our 
13 CO and C I8 maps prevent any robust analysis without addi- 
tional shorter spacing data. 



3. Results 

3.1. SMA 1.3-mm dust continuum map 

The 1.3-mm continuum map shown in the left hand panel of 
Figure [T] shows a total of three sources detected above the 10-<x 
level. The two central sources are coincident with the two com- 
ponents (VLA 1623 to the east, Knot-A to the west) previously 
detected at both 2. 7-mm ( Looney et al. 2000} and 7-mm ( Ward- 
Thompso n et al.|[201 l| l. The third, westernmost source is coinci- 
dent with the cm source Knot-B by Bontemps & Andre ( 1997}. 
The separation between the two central continuum sources is 
1.2", i.e -1 70 AU assuming a distance D=139±6 pc for the p 
Oph doug* 3 ] ( Mamaj ek| |2008j ) . Knot-B is located 10.4" away of 
VLA1623, i.e -1450 AU toward the west. 

Table 1 gives the positions and deconvolved sizes (FWHMs) 
of the millimeter continuum sources in our SMA map, obtained 
by fitting the continuum visibilities (in the uv-plane) with a 
model that includes emission from three 2-D Gaussian sources 
and shown in Fig. [2] This modeling work confirms that the pro- 
file of VLA 1623 is typical of an envelope profile, suggesting that 
this component dominates the single-dish large-scale continuum 
emission and is representing the centre of the protostellar enve- 
lope of VLA1623. The continuum profiles of Knot-A and Knot- 
B, shown in Fig. [2] show almost no dependency with baseline 
length, which is typical of much more compact sources. 

The three sources have measured peak and integrated fluxes 
as reported in Table 1. The total integrated 1.3-mm emis- 
sion above the 3<x level on the whole primary beam field is 
300±80 mjy. We recover 20% of the total integrated flux of 
VLA1623 as seen with t he IRAM-30m (15 00 mjy within a re- 
gion of diameter 20", see Andre et al. 1993[ ). Such a low percent- 
age is not surprising considering that our high-resolution obser- 
vations do not probe the large spatial scales, and therefore filter 
out most of the extended envelope emission. 



3.2. SMA l2 CO(2-1) map 

Since the shortest baseline during our SMA observations was 
-68 m (corresponding to a maximum spatial scale -4.7" probed 



3 VLBA (Loinard et al. 2008) and extinction jLombardi etaT 2008 1 
methods for distance determination are in favor of a slightly smaller 
distance (~120 pc) for the p Oph core. 
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Fig. 1. Left panel: Wide-field SMA 1 .3-mm continuum map towards VLA1623. In this map the beam size is 0.7"x0.4", allowing the 
three components detected at the positions of VLA1623, Knot-A and Knot-B to be disentangled. The rms noise in the continuum 
map is 2.1 mJy/beam. Black contours are levels of 6.3 mJy/beam (3cr-level, dashed contour), then 5<x to 20<x levels by 5cr, and 
ultimately 30<r and 40<x levels. The two white arrows represent the opening angle (~ 30°) of the protostellar outflow according to 
CO low-resolution observations and knots of shocked H2 emission ( [Andre et al. 1990 Dent et al. 1995| l. - Right Panel: Blow-ups 
of the 1.3-mm continuum map towards the three detected sources. The positions (blue stars) and FWHMs (dashed ellipses) of the 
deconvolved sources modeled (see Table 1) are shown. 
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Fig. 2. The left panel shows, for the three continuum sources, the best-fit model profiles of the visibilities' real parts as a function of 
projected baseline length. The visibilities for the continuum source modeled at the position of VLA1623 are shown as green dots, 
the model visibilities towards Knot-A as blue crosses, and the model visibilities reproducing the source at the position of Knot-B 
are shown as open orange squares. The right panel shows the observed profile of the visibilities real parts in our SMA continuum 
data (plotted as black points with 3cr error bars) centered at the position of VLA1623, and the global best-fit model (red crosses) 
featuring the three gaussian sources shown in the left panel, centered at the same position. The properties of the three continuum 
sources deduced from this best-fit model are reported in Table 1 . 



in the image plane), our 12 CO(2-l) map is affected by interfer- 
ometric spatial filtering. For example, our SMA map recovers 
only ~1% of the high-velocity CO(2-l) emission mapped by 
Andre et aL| ( [T9 90) with their single-dish observations. Thus, not 
surprisingly, we do not detect the large-scale protostellar outflow 
observed by |Andre et al. ( 1990[ ), nor the extended emission from 
the envelope, which are almost entirely resolved out. We also 
suspect that interferometric filtering and/or self-absorption, sim- 
ilar to what is seen in the single-dish observations of VLA1623, 
are responsible for the lack of CO emissi on in our map near 
the systemic velocity (vlsr ~3.8 kms' 1 , |Narayanan & Logan 
2006). Therefore, we stress that the structure seen in our CO map 
can only trace emission from compact material, at redshifted and 
blueshifted velocities. 

While tentative (~3<x detection) 12 CO emission is detected 
towards VLA1623 (see Fig(3^), at slightly blueshifted velocities 
(~2kms _1 ), the only prominent (detections above the 5<x level) 
emission features in our CO map are two small-scale sources, 
closely associated with Knots A and B. Figures [3]; and|4i show 
the CO integrated emission contours towards Knot-A and Knot- 
B. In both cases, the CO emission is located right beyond the 



continuum peak, ~0.2-0.4" downstream along the propagation 
direction of the VLA1623 outflow. The integrated spectra shown 
in Figs[3jl and|4]3 reveal that, while the CO emission associated 
with Knot-A has two components (red- and blueshifted) at low 
projected velocities with respect to the systemic velocity of the 
cloud, only blueshifted CO emission is associated with Knot-B. 

Our results therefore suggest that, in a region ~20"x20" 
around VLA1623, the brightest compact 12 CO emission struc- 
tures (shifted from the systemic velocity of the cloud) are asso- 
ciated with Knot-A and Knot-B. 



4. Nature of the detected sources 

Our SMA data suggest that, while the 1.3-mm continuum 
emission from Knot-A and Knot-B is almost point-like, the 
continuum emission from VLA1623 is partially extended and 
its visibility curve resembles envelope profiles seen on larger 
scales (|Andre et al.||199"3) |Ohashi et al.||199"71 see electronic 
Figure|2]). This, and the fact that the continuum emission towards 
VLA1623 presents the greatest integrated flux, offer further ev- 
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Fig. 3. 12 CO(2-l) emission towards the center of the VLA1623 
envelope. Significant CO (above the 10cr level) emission, asso- 
ciated with the continuum source Knot- A, is detected, (a) and 
(b) show the peak spectrum towards VLA1623 and the peak 
spectrum of the emission feature associated with Knot-A, re- 
spectively. In all panels, the central black dotted vertical line 
indicates the systemic velocity v l $r ~3.8 kms -1 , while the col- 
ored dashed vertical lines show the positions of the redshifted 
and blueshifted components seen in the 12 CO(2-l) single- 



dish spectrum towards VLA1623 ( | Andre et al. 1990 Dent 
et al. 1995] >. (c) shows the 1.3-mm continuum emission map, 
and contours of 12 CO(2-l) integrated emission: blueshifted 
from -6kms~' to +3.0kms _I , redshifted from +4.5kms~ 1 to 
+ 10 km s~'). The rms noise in the integrated CO emission map is 
cr ~0.15 Jy/beamkm s _1 , and the contours shown are levels of 3- 
cr to 1 1-cr by 3-0" (~3cr dashed contours). The black crosses and 
the white squares mark the positions of the continuum sources 
at 2.7-mm (Looney et al. 2000[ ) and 7-mm ( [Ward- T hompson 
et al. 201 \) , respectively, (d) shows the CO spectrum, integrated 
over a 3" box, including the whole extent of continuum emission 
from both VLA1623 and Knot-A. 



idence that VLA1623 traces the center of the protostellar enve- 
lope seen on larger scales. 

The close alignment of both Knot-A and Knot-B with the 
southwestern outflow cavity wall, as traced by several previous 
studie s ( [Andre et al.| [19901 |Dent et al.| [1993) |Yu & Chernin[ 
|1997| l, along with their compactness at millimeter wavelengths 
and their association with centimeter free-free emission first sug- 
gests that these two sources are outflow features, i.e shocks be- 
tween the outflowing material from VLA1623 and the dense am- 
bient material, rather than protostellar components. Indeed, it 
has been suggested (Lee et al. 2001 Cunningham et al. 2005) 
that the interaction between circumstellar material and outflow- 
ing material is likely to create high-density regions in the cavity 
walls, where the development of a shear layer could also be the 
site of maser and free-free emission. Interpretation of Knot-A 
as a protostellar companion of VLA1623, which is based on the 
detection of 2.7-mm continuum emission towards both sources 
by |Looney et al.| ( [20 00), has been proposed when no millimeter 
continuum emission towards Knot-B had been reported yet. We 
stress that, while it would still be statistically possible to detect 
one protostellar companion (such as Knot-A) aligned with the 
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Fig. 4. 12 CO(2-l) emission towards Knot-B. Significant 12 CO 
(above the ~5cr level) emission is detected, as a compact source 
~0.2" away from the 1.3-mm continuum source tracing Knot- 
B. (a) shows the 1.3-mm continuum emission map, and the 
I2 CO(2-l) integrated emission contours (blue from -8kms~' 
to +3.0 km s -1 , red from +4.5kms _I to +10kms~ 1 ). The solid 
contours are levels of 3<x and 5tr, and the dashed contours show 
the levels of -3<x, the rms noise being of 0.19 Jy/beamkm s _1 . 
(b) shows the integrated 12 CO spectrum of the emission feature 
associated with Knot-B. The middle black dashed vertical line 
shows the systemic velocity Vlsr ~3.8 kms -1 , while the col- 
ored dashed vertical lines show the position of the redshifted and 
blueshifted components seen in the 12 CO(2-l) single-dish spec- 
trum towards VLA1623 ( [Andre et al.| \1§9U[ |Dent et al.| \T§§5\ . 
The noise level is greater towards the spectrum of Knot-B than 
towards VLA1623 and Knot-A, because this source is located 
~10" away from the phase center chosen for our observations. 



cavity wall of the main protostellar source, the detection of two 
perfectly aligned companions along this axis (such as Knot-A 
and Knot-B) is highly unlikely. 

Both Knot-A and Knot-B are associated with small-scale, 
low-velocity CO emission in our SMA map, detected at a po- 
sition shifted by ~0.2" to ~0.4" from the continuum emission 
(downward along the outflow propagation direction). Such sys- 
tematic shift between the 1.3-mm continuum and CO emission 
peaks does not favor a scenario where the CO emission would 
come from a small disk surrounding a young stellar object, be- 
cause such a disk would also be traced by dust continuum emis- 
sion at the position of the CO peak. 

While the CO single-dish observations of VLA1623 trace 
only redshifted velocities for the gas towards the western outflow 
lobe on large scales, the two CO sources associated with Knot- 
A and Knot-B mostly show blueshifted CO emission. Murillo & 



Lay (2010[) also report detecting blueshifted CO(3-2) emission 



along the same redshifted southwestern outflow cavity wall from 
VLA1623. Some detections of reverse velocity gas along out- 
flow cavity walls have been reported (see for example the case 
of Serp-SMMl by van Kempen et al.|20 09b ). These features can 
be explained as a result of the development of Kelvin-Helmholtz 
and Rayleigh-Taylor instabilities in outflow cavities (see for ex- 
ample the hydrodynamical simulations by [Cunningham et al. 
2005 ). Another explanation would be that the outflow driven by 
VLA1623 has an axis close to the plane of the sky, as already 
proposed by Andre et al. ( 1990 1 to explain its low projected ve- 
locities, and that it opens a cavity that has a larger opening angle 
than its angle with the plane of sky. Such a configuration would 
lead to overlapping redshifted and blueshifted CO emission de- 
tected on both sides of the source, due to the projection of the 
front and back edges of the cavity lobe along the line of sight 
dCabrit & Bertout[[T986l l. 

The SMA CO sources that we detect are elongated in the per- 
pendicular direction with respect to the outflow propagation di- 
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rection. This is especially visible in the case of the CO emission 
associated to Knot-A (FWHM 1. 1"x0.5", P.A 35°), because it is 
more extended. Moreover, the CO emission associated to Knot- 
B has a curved shape, with an apex oriented towards the jet's 
propagation direction (see Figure [4]}), resembling a bow-shock 
shape. This suggests that the 12 CO compact emission detected 
towards Knots A and B does not trace outflowing material, but is 
instead due to material interacting with the outflow. Such spatial 
properties are consistent with a picture where the 12 CO sources 
originate in small-scale C-type shocks along the walls, as de- 
scribed in the model by Visser et al. ( 201 1\ and already sug- 
gested by van Kempen et al. (2010) towards the outflow cavity 
ofHH46. 

Furthermore, the brightness temperature of these CO small- 
scale sources (~60 K towards Knot-A, 30 K towards Knot-B, see 
Figs.[3]and|4]i is higher than the CO brightness temperature from 
the large-scale outflowing material probed by single-dish obser- 
vations (~20 K). Such an increase in the CO brightness temper- 
ature in outflow cavities is generally interpreted as a signpost of 
gas heating, due either to outflow shocks (Bachiller & Tafalla 



1999; Hollenbach 1997) or to ultraviolet photons produced by 



the protostellar accretion ( |van Kempen e t al. 2009a) 

Following the new compelling arguments brought by our 
SMA observations of VLA1623 presented above, we propose 
that Knot-A and Knot-B might represent outflow-features, due 
to shocks between the outflowing gas and the ambient material 
along the southwestern cavity wall. Emission from outflow fea- 
tures in the millimeter continuum maps of young protostars has 
already been reported by lGueth et al. (|2003 1 towards LI 157, Yen 



et al. ( 2010)> towards B335, or|van Kempen et al.|(|2009a|> towards 



HH 46, on larger scales. While further observations of these 
sources are needed to assess the exact nature of the millime- 
ter continuum emission, it is likely that the emission is due to a 
blend of free-free emission (from the shock) and dust continuum 
emission (from the local accumulation of material). The bright 
CO emission detected in Knots A and B might also come from 
column density enhancement and/or heating of the gas, similar 
to what was observed towards IRAS 16293 by|Yeh et al.| p008), 
and towards HH 46 by van Kempe net al.| ( |2010| l thanks to high-J 
CO lines emission. Complementary molecular observations are 
needed to confirm the exact physical processes (shocks, photon 
heating) responsible for both the CO and continuum emission 
towards these sources. 

These observations provide the clearest example so far of 
outflow/envelope interactions in the inner few hundred AU of 
protostars, resulting in compact dust continuum emission. Given 
the importance of outflows and jets during the early protostellar 
stages, such features may be common. It therefore demonstrates 
the need for thorough analysis of the millimeter continuum com- 
ponents detected close to the axes of protostellar jets, since they 
might be due to outflow features rather than to protostellar com- 
panions. 

Finally, while it was considered as a good example of a 
close protobinary system so far, the present analysis suggests 
that VLA1623 is single on the scales a > 100 AU probed by 
our SMA observations. This supports the conclusion of Maury 
|et al.| ( [20T0| , who propose a low rate of small-scale (1000 AU> 
a >100 AU) multiplicity for Class protostars, also because 
VLA1623 is the prototype for protostars at the Class stage. 
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